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Abstract 



Study of fine spectral properties of quasiperiodic and similar discrete 
Schrodinger operators involves dealing with problems caused by small de- 
nominators, and until recently was only possible using perturbative methods, 
requiring certain small parameters and complicated K AM- type schemes. We 
review the recently developed nonperturbative methods for such study which 
lead to stronger results and are significantly simpler. Numerous applications 
mainly due to J. Bourgain, M. Goldstein, W. Schlag, and the author are also 
discussed. 
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and Vin.m) = VnS{n,m) is a potential given by K = f{T^' ■ ■■T^-'O), € T\ 



where Ti9 = 9 + LUi, and co is an incommensurate vector. In certain cases A may also 
be replaced by a long-range Laplacian. Replacing T^'s with other commuting ergodic 
transformations would give a general framework of ergodic Schrodinger operators 
([21]; see Sec. 8. for an example of this kind) but we will mostly focus on the 
quasiperiodic (QP) operators that have been intensively studied in Physics and 
Mathematics literature. For another review of some recent developments in this 
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1. 



Introduction 

Consider an operator acting on on £^(Z'') defined by 



Hx^A + XV, 



(1.1) 



where A is the lattice tight-binding Laplacian 
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area sec [7] . The questions of interest are the nature and structure of the spectrum, 
behaviour of the eigenfunctions, and particularly the quantum dynamics: properties 
of the time evolution = e'*^^'o of an initially localized wave packet ^q. 

Of particular importance is the phenomenon of Anderson localization (AL) 
which is usually referred to the property of having pure point spectrum with ex- 
ponentially decaying eigenfunctions. A somewhat stronger property of dynamical 
localization (see Sec. 7.) indicates the insulator behavior, while ballistic transport, 
which for = 1 follows from the absolutely continuous (ac) spectrum, indicates the 
metallic behavior. 

Operators with ergodic potentials always have spectra (and pure point (p.p.) 
spectra, understood as closures of the set of eigenvalues) constant for a.e. realization 

of the potential. Moreover, the p.p. spectrum of operators with ergodic potentials 
never contains isolated eigenvalues, so p.p. spectrum in such models is dense in a 
certain closed set. An easy example of an operator with dense pure point spectrum 
is which is operator (1.1) with A^^ = 0, or pure diagonal. It has a complete 
set of eigenfunctions, characteristic functions of lattice points, with eigenvalues Vj. 
H\ may be viewed as a perturbation of i?oo for small A~^. However, since Vj are 
dense, small denominators (Vi — Vj)^^ make any perturbation theory difficult, e.g. 
requiring intricate KAM-type schemes. 

The probabilistic KAM-type scheme was developed by Frohlich and Spencer 
[26] for random potentials {Vn are i.i.d.r.v.'s) in the multi-dimensional case, and 
is called multi-scale analysis. It was significantly modified, improved, and widely 
applied in the later years by a number of authors, most notably [24] . An alternative 
method for random localization was found by Aizenman-Molchanov [2] and later 
further developed by Aizenman and coauthors. While still requiring certain large 
parameters this method relies on direct estimates of the Green's function rather 
than a step-by-step perturbation scheme. 

For QP potentials none of the above methods work, as, among other reasons, 
they do not allow rank-one perturbations, nor Wegner-type estimates. The situation 
here is more difficult and the theory is far less developed than for the random case. 
With a few exceptions the results are confined to the ID case, and also 1-frequency 
case (6=1) has been much better developed than that of higher frequencies. 

One might expect that H\ with A small can be treated as a perturbation 
of Hq = A, and therefore have ac spectrum. It is not the case though for random 
potentials in d — 1, where AL holds for all A. Same is expected for random potentials 
in d = 2 (but not higher). Moreover, in ID case there is strong evidence (numerical, 
analytical, as well as rigorous [8]) that even models with very mild stochasticity in 
the underlying dynamics have point spectrum for all values of A like in the random 
case (e.g. Vn = Xf{n'^a + 6) for any a > 1). At the same time, for QP potentials 
one can in many cases show ac spectrum for A small as well as pure point spectrum 
for A large (sec below), and therefore there is a metal insulator transition in the 
coupling constant. It is an interesting question whether quasiperiodic potentials are 
the only ones with metal-insulator transition in ID. 
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2. Perturbative vs nonperturbative 

It is probably fair to say that much of the theory of qusiperiodic operators has 
been first developed around the almost Mathieu operator, which is 

Hx,u.,e = ^ + Xf{e + nw) (2.1) 

acting on ^^(Z), with / : T ^ T; f{9) = cos{2'k0). The first KAM-type approaches, 
in both large and small coupling regimes, were developed for this or similar models 
[23, 5]. The perturbative proofs of complete Anderson localization for A > A(w) 
large are due to Sinai [41] and Prohlich-Spencer-Wittwer [27], and both applied 
to cos-type /. For A small Chulaevsky-Delyon [19] proved pure ac spectrum using 
duality and the construction of Sinai [41]. Elliasson (see [25] for a review) developed 
alternative KAM-type arguments for both large and small A for the case of real- 
analjrtic (actually, somewhat more general) class of / in (2.1). 

The common feature of the perturbative approaches above is that, besides 
all of them being rather intricate multi-step procedures, they rely extensively on 
eigenvalue and eigenfunction parametrization and perturbation arguments. 

The common feature of the perturbative results in the quasiperiodic setting 
is that they provide no explicit estimates on how large (or small) the parameter 
A should be, and, more importantly, A clearly depends on w at least through the 
constants in the Diophantine characterization of lo. 

In contrast, the nonperturbative results allow effective (in many cases even 
optimal) and, most importantly, independent of uj, estimates on A. We will take the 
latter property (uniform in lo estimates on A) as a definition of a NP result. 

Recently developed nonperturbative methods are also quite different from the 
perturbative ones, in that they do not employ multi-scale schemes: usually only 
a few (from one to three) sufficiently large scales are involved, do not use the 
eigenvalue parametrization, and rely instead on direct estimates of the Green's 
function. They are also significantly less involved, technically. One may think that 
in these latter respects they resemble the Aizenman-Molchanov method for random 
localization. It is, however, a superficial similarity, as, on the technical side, they 
are still closer to and do borrow certain ideas from [26, 24] . 

Several results that satisfy our definition of nonperturbative appeared prior 
to the recent developments, and were all related to the almost Mathieu operator 
(see [32] for a review). In [33, 34] AL was proved for A > 15, and existence of p.p. 
component for A > 2. The latter papers, while introducing some of the ingredients 
of the recent nonperturbative methods, did not take advantage of the positivity of 
the Lyapunov exponents which proved very important later. 

3. Lyapunov exponents 

Here for simplicity we consider the quasiperiodic case, although the definition 
of the Lyapunov exponents and some of the mentioned facts apply more generally 
to the ID ergodic case. 
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For an energy E gM. the Lyapunov exponent 7(-B) is defined as 



lim 



J^\n\\Mk{e,E)\\de 



(3.1) 



n- 



k 



where 





Mk{e,E)= JJ 



'E-Xf{ujn + e) -1 
1 



) 



n—k—1 



is the fc— step transfer-matrix for the eigenvalue equation H'^ = E^. 

We will be interested in the regime when Lyapunov exponents are positive 

for all energies in a certain interval intersecting the spectrum. It is well known 
and fairly easy to see that if this condition holds for all E e 9., there is no ac 
component in the spectrum for a.e. 6 (it is actually true for all 9 [39]). Positivity of 
Lyapunov exponents, however, does not imply exponential decay of eigenfunctions 
(in particular, not for the Liouvillc lo [3] nor for the resonant 9 gT'' [36]). 

NP methods, at least in their original form, stem to a large extent from esti- 
mates involving the Lyapunov exponents and exploiting their positivity. 

The general theme of the results on positivity of j(E), as suggested by pertur- 
bation arguments, is that the Lyapunov exponents are positive for large A. This was 
first established by Aubry-Andre [1] for the almost Mathieu operator with A > 2. 
Their proof was made rigorous in [3]. Another proof, exploiting the subharmonic- 
ity, was given by Herman [31], and applied to trigonometric polynomials /. The 
lower bound in [31] was in terms of the highest coefficient of the trigonometric 
polynomial and therefore this did not easily extend to the real analytic case. All 
the subsequent proofs, however, were also based on subharmonicity. Sorets-Spencer 
[42] proved that for nonconstant real analytic potentials ?; on T (6 = 1) one has 
jiE) > ^InA for A > X{v) and all irrational oj. Another proof was given in [11], 
where this was also extended to the multi- frequency case (6 > 1) with, however, 
the estimate on A dependent on the Diophantine condition on oj. Finally, Bourgain 
[10] proved that Lyapunov exponents are contimious in lu at every incommensurate 
u) (for 6 > 1; for 6 = 1 this was previously established in [16]), and that led to the 
following Theorem which is the strongest result in this general context up to date: 

Theorem 1 [10] Let f be a nonconstant real analytic function on T^, and H given 
by (1.1). Then, for A > A(/), we have ^{E) > 5 In A for all E and all incommen- 
surate vectors u). 

3.1. Corollaries of positive Lyapunov exponents 

The almost Mathieu operator On one hand the almost Mathieu operator, 
while simple-looking, seems to represent most of the nontrivial properties expected 
to be encountered in the more general case. On the other hand it has a very special 
feature: the duality (essentially a Fourier) transform maps Hx to H^/Xi hence A = 2 
is the self-dual point. Aubry-Andre [1] conjectured that for this model, for irrational 
OJ a sharp metal-insulator transition in the coupling constant A occurs at the critical 
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value of coupling A = 2 : the spectrum is pure point for A > 2 and pure ac for A < 2. 
A second, related, conjecture was that the dual of ac spectrum is pure point and 
vice versa. Both conjectures were modified based on the results of [3, 37, 36]. The 
first modified conjecture stated pure point spectrum for Diophantinc uj and a.e. 
for A > 2 and pure ac spectrum for A < 2 for all lo,9. As for the duality, the 
question, after some prior developments, was resolved in [28] where it was shown 
that the dual of point spectrum is ac spectrum (the proof applied in a more general 
context), and it was used (together with [38, 30]) to prove that the spectrum is 
purely singular continuous at A = 2 for a.e. w, 9. 

As with the KAM methods, the almost Mathicu operator was the first model 
where the positivity of Lyapunov exponents was effectively exploited: 

Theorem 2 [32] Suppose oj is Diophantine and ^{E^u)) > for all E G [_Ei,i?2]. 
Then the almost Mathieu operator has Anderson localization in [Ei,E2\ for a.e. 9. 

The condition on 6 in [32] was actually explicit (arithmetic) and close to op- 
timal. This, combined with the mentioned results on the Lyapunov exponents for 
the almost Mathieu operator [31] and duality [28] led to the following corollary: 

Corollary 3 The almost Mathieu operator H^^x^e has 

1° [32] for A > 2, Diophantine u and almost every 9 gM., only pure point 

spectrum with exponentially decaying eigenf unctions. 
2" [28] for A = 2, and a.e. uj,9 gM. purely singular- continuous spectrum. 
3° [32, 28] for A < 2, Diophantine w G R and a.e. 9 £ M, purely ac spectrum. 

Precise arithmetic descriptions of u), 9 are available. Thus the Aubry- Andre 
conjecture is settled at least for almost all oj, 9. One should mention, however, that 
while 1° is almost optimal, both 2° and 3° are expected to hold for all 6 and all 
w ^ Q, and such extension remains a challenging problem (see [40]). 

The method in [32], while so far the only nonperturbative available allowing 
precise arithmetic conditions, uses some specific properties of the cosine. It extends 
to certain other but rather limited situations. A much more robust method was 
developed by Bourgain-Goldstein [11], which allowed them to extend (a measure- 
theoretic version of) the above result to the general real analytic as well as the 
multi-frequency case. Note, that essentially no results were previously available for 
the multifrequency case, even perturbative. 

Theorem 4 [11] Let f be non-constant real analytic on T*" and H given by (2.1). 
Suppose ^{E,oj) > for all E s [EijE^] and a.e. w € T^. Then for any 9, H has 
Anderson localization in [Ei,E2] for a.e. w. 

Combining this with Theorem 1 one obtains [10] that for A > A(/), H as above 
satisfies Anderson localization for a.e. w. 

One very important ingredient of the method of [11] is the theory of semi- 
algebraic sets that allows one to obtain polynomial algebraic complexity bounds for 
certain "exceptional" sets. Combined with measure estimates coming from the large 
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deviation analysis of ^ In ||M„(0)|| (using subharmonic function theory and involv- 
ing approximate Lyapunov exponents), this theory provides necessary information 
on the geometric structure of those exceptional sets. Such algebraic complexity 
bounds also exist for the almost Mathicu operator [32] and are actually sharp albeit 
trivial in this case due to the specific nature of the cosine. 

4. Without Lyapunov exponents 

While having led to significant advances, Lyapunov exponents have obvious 
limitations, as any method, based on them, is restricted to ID nearest neighbor 
Laplacians. It turns out that the above methods can be extended to obtain NP 
results in certain quasi-lD situations where Lyapunov exponents do not exist. 

For the next Theorem let H be an operator (1.1) defined on £"^{8) where 
S = 'ZiX So, is a strip. 5*0 here is a finite set with a metric, and dist((n, s), {n' , s')) = 
\n - n'\ + dist(s, s'). Let V(n,s) = fsiO + nw), 6* e T. 

Theorem 5 [14] Assume fs, s G 5*0, are non-constant real analytic functions on 
T. Then for any 6 & T and A > A(/s), operator H has AL for a.e. co. 

The following nonperturbative Theorem deals with the case of small coupling: 

Theorem 6 [15] Let H be an operator (2.1), where f is real analytic on T and uj 
is Diophantine. Then, for A < A(/), H has purely ac spectrum for a.e. 0. 

We note that an analogue of this Theorem does not hold in the multi- frequency 
case (see next section). Theorem 6 is a result on non-perturbative localization in 
disguise as it was obtained using duality [28] from a localization Theorem for a 
dual model which has in general a long-range Laplacian and was in turn obtained 
by an extension of the method of [32]. A certain measure-theoretic version of it 
by the method of [11] allowing non-local Laplacians but leading only to continuous 
spectrum is also available [6]. Theorem 5 was obtained by an extension of the 
method of [11]. Both Theorems above rely on large deviations for the quantities 
of the form iln|det(-ff — E)\\ and path-determinant expansion for the matrix 
elements of the resolvent [15]. The methods developed in [15] apply also to certain 
other situations with long-range Laplacians, for example the kicked rotor model 
(see Sec. 8.). 

5. Multidimensional case: d > 1 

As mentined above, there are very few results in the multidimensional lattice 
case {d> 1). Essentially, the only result that existed before the new developments 
was a perturbative Theorem - an extension by Chulaevsky-Dinaburg [20] of Sinai's 

[41] method to the case of operator (1.1) on I'^^lf) with Vn = A/(n • w), w S R'^, 
where / is a cos-type function on T. Recently, Bourgain [6] obtained this result for 
real analytic / by a nonperturbative method. Note that since 6=1, this avoids 
most serious difficulties and is therefore significantly simpler than the general multi- 
dimensional case. 
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Theorem 7 [20, 6] For any e > there is A(/, e), and, for \ > A(/, e), n{X, /) C T 
with mes{fl) < e, so that for co ^ fl, operator (1.1) with Vn as above has Anderson 
localization. 

This should be confronted with the following Theorem of Bourgain [8] 

Theorem 8 [8] Let d = 2 and f{9) = cos27r6' in H = defined as above. Then 
for any A measure of oj s.t. H^j has some continuous spectrum is positive. 

Therefore for large A there will be both iv with complete^ localization as well 
as those with at least some continuous spectrum. This shows that nonperturbative 
results do not hold in general in the multi-dimensional case! 

A similar (in fact, dual) situation is observed for ID multi- frequency {d = 
1; 6 > 1) case at small disorder. One has, by duality: 

Theorem 9 [20, 19] Let H be given by (2.1) with 0, w S T'' and f real analytic on 
. Then for any e > there is A(/, e) s.t. for A < A(/, e) there is fl{X, f) C T*" 
with mes{^l) < e so that for to ^ fl, H has purely ac spectrum. 

Theorem 10 [8] Let d = 1, 6 = 2 and f be a trigonometric polynomial on with 
a non-degenerate maximum. Then for any A measure of oj s.t. has some point 
spectrum, dense in a set of positive measure, is positive. 

Therefore, unlike the 6=1 case (see Theorem 6), nonperturbative results do 
not hold for absolutely continuous spectrum at small disorder. 

6. Perturbative results by NP methods 

While the above demonstrates the limitations of the NP results, the nonper- 
turbative methods have been applied to significantly simplify the proofs and obtain 
new perturbative results that previously have been completely beyond reach. 

We refer the reader for a description of many such applications to [7, 6] . In par- 
ticular, new results on the construction of QP solutions in Melnikov problems and 
nonlinear PDE's, obtained by using certain ideas developed for NP quasi-periodic 
localization (e.g. the theory of semi- algebraic sets) are presented there. 

We will only mention here a theorem by Bourgain-Goldstcin-Schlag that is 
the only one so far treating a "true" d > 1 situation. Note that here d = 2, and 
the reasons why it has not yet been extended to higher dimensions are not just 
technical, but conceptual (there are certain purely arithmetic difficulties). 

Theorem 11 [12] Let d = b = 2 and f be real analytic on T"^ such that all functions 
/(^i) ■)> /(■> ^2), (^1, ^2) € are nonconstant. Then for any e > there is \{f, e) 
s.t. for A > A(/, e) there is J7(A, /) C T'' with mes{n) < e so that for u> ^ Q, 
operator (1.1) with Vn = A/(nia;i, n2W2) has Anderson localization. 
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7. Dynamical localization 

Anderson localization does not in itself guarantee absense of quantum trans- 
port, or nonspread of an initially localized wave packet, as characterised, e.g., by 
boundedness in time of moments of the position operator [22] ([35] for an example 
of physical model with coexistence of exponential localization and quantum trans- 
port). Considering for simplicity the second moment 

^ Jo n 

we will say that H exhibits dynamical localization (DL) if {x^)t < Const. We will 
say that the family {Hq} exhibits strong DL if Jr^b dO a\xp{x^)t < Const. We 

note that the results mentioned below will hold with more restrictive definitions of 
DL (involving the higher moments of the position operator) as well. DL implies 
p.p. spectrum by RAGE theorem (see, e.g. [21]), so it is a strictly stronger notion. 

It turns out that nonperturbative methods allow for such dynamical upgrades 
as well. For the almost Mathieu operator we have 

Theorem 12 [29] For A > 2 and Diophantine u (as in [32]) strong DL holds. 

While proved in [29] with a slightly more restrictive condition on cj. Theorem 
12 holds as stated by a result of [16]. For the results obtained by methods stemming 
from the approach of [11] one has 

Theorem 13 [14] In Theorems 4,5,7,11 dynamical localization also holds. 

This also applies to other results on nonperturbative localization, e.g. [13]. 



8. Quantum kicked rotor 

The quantum kicked rotor was introduced in [18] as a model in quantum chaos. 
It is given by the time-dependent Schrodinger equation on I/^(T) 

where y(t, &) = KCOs27r^^^^^ 5{t — n). It represents quantization of the Chirikov 
standard map, and a conjecture (e.g. [4]) was that for a.e. a,b the solution is 
almost-periodic in time, thus demonstrating "quantum suppression of chaos" (as 
some chaos is expected for the standard map) . Such almost periodicity follows from 
dynamical localization for the Floquet operator W : L^(T) L^(T) defined by 
W4'{t, 6) — (p{t + 1, 6'). H'^ is a unitary operator that in Fourier representation can 
be written as a product U ■ S where S{n, m) = S{n — m) is a Toeplitz operator 
with very fast decay of S{n) and D{n,m) = D{n)S{n,m) is a diagonal operator 
with D{n) = cxp(27ri(T".x)2), with T being the skew shift of the torus T{xi,X2) = 
{xi + (jj,X2 + xi), and xi,X2,(jJ determined by a,b. The nonperturbative methods 
(particularly, the method for skew-shift dynamics localization in [13] and a long- 
range method [15] ) were further developed for this model to obtain 
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Theorem 14 [9] For any e > and any fixed b there is K{e), and for k < K{e), 
Q{k) C T with mes{fl) < e, s.t. for a ^ ^{k), operator W has DL. 

Exploiting the multiplicative nature of W one also obtains a nonperturbative 

counterpart: 

Theorem 15 [17] There is kq > such that for any b, operator W has dynamical 
localization for k < kq and a.e. a. 

This confirms the "quantum suppression of chaos" conjecture for small k. 
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